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ABSTRACT 

In the development and testing of novel structural and controls concepts, such as morphing aircraft wings, appropriate 
models are needed for proper system characterization. In most instances, available system models do not provide the 
required additional degrees of freedom for morphing structures but may be modified to some extent to achieve a 
compatible system. The objective of this study is to apply wind tunnel data collected for an Unmanned Air Vehicle 
(UAV), that implements trailing edge morphing, to create a non-linear dynamics simulator, using well defined rigid 
body equations of motion, where the aircraft stability derivatives change with control deflection. An analysis of this 
wind tunnel data, using data extraction algorithms, was performed to determine the reference aerodynamic force and 
moment coefficients for the aircraft. Further, non-linear influence functions were obtained for each of the aircraft’s 
control surfaces, including the sixteen trailing edge flap segments. These non-linear controls influence functions are 
applied to the aircraft dynamics to produce deflection-dependent aircraft stability derivatives in a non-linear dynamics 
simulator. Time domain analysis of the aircraft motion, trajectory, and state histories can be performed using these non- 
linear dynamics and may be visualized using a 3-dimensional aircraft model. Linear system models can be extracted to 
facilitate frequency domain analysis of the system and for control law development. The results of this study are useful 
in similar projects where trailing edge morphing is employed and will be instrumental in the University of Maryland’s 
continuing study of active wing load control. 
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1. INTRODUCTION 

In the development and testing of novel structural and controls concepts, such as morphing aircraft wings, appropriate 
models are needed for proper system characterization. These models provide a means for understanding the 
characteristics and behaviors of such concepts in real world applications, but may be difficult to derive. In this study, a 
well understood aircraft dynamics model was modified to represent an aircraft with non-conventional control surfaces. 

Development of an aircraft dynamics simulation requires the gathering of information about the aircraft in question. 
Extensive wind tunnel data was used to obtain the aerodynamics force, moment, and damping derivatives for an 
Unmanned Aerial Vehicle (UAV) with main wing trailing edge segmented flaps. This data was used to create a 
component build up method, in the form of non-linear controls influence functions, for determining an aircraft’s 
stability derivatives at any point in the simulation, allowing for stability derivatives that change with time and control 
deflection magnitude. These stability derivative models are implemented in a non-linear aircraft dynamics simulation. 

The motivation for this paper comes from the need to understand the dynamics of aircraft with non-conventional control 
surfaces, such as morphing aircraft. These aircraft provide a wide range of new capabilities, in the form of additional 
degrees of freedom, which can be used to perform functions in addition to providing conventional control. The 
resulting aircraft dynamics simulation provides a multiple functionality tool that can be used for aircraft performance 
characterization and the extraction of system models for controls design. This study is in support of the University of 
Maryland’s Lift Distribution Control (LDC) research. 



2. BACKGROUND AND METHODOLOGY 


2.1. Aircraft equations of motion 

The aircraft equations of motion used to build this model are well-defined equations. Stevens and Lewis 1 , and Pamadi 2 
provide two sources for in-depth understanding of aircraft dynamics and are left as references to the conventional 
aircraft dynamics notation used in this paper. The rigid-body non-linear equations of motion, as presented in Stevens 
and Lewis 1 , and referenced to the aircraft body-fixed frame, are given by equations (1) - (12). 


Force equations: 
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Linearized models are extracted from these non-linear equations and their standard form must be used. The equations 
of motion are rearranged and rotated into other coordinate frames, such as the wind coordinate frame, to organize the 
aircraft state vector into the standard aerospace form: 
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2.2. Testbed aircraft 

This study is concerned with understanding the dynamics of a UAV with novel control effectors. The aircraft, a former 
US Army target drone classified as the FQM-117B, extensively modified, housed, and flown by NASA Langley 
Research Center (LaRC), implements segmented trailing edge flaps to achieve trailing edge morphing. Segmented 
trailing edge flaps consist of many independently controlled aileron segments that provide additional degrees of 
freedom in aircraft control. This one-of-a-kind aircraft is part of NASA LaRC’s FLying Controls testbed (FLiC) 
program used for neural network and other advanced controller development and testing. Figure 1 shows these 
segmented flaps and their corresponding identification number. 




Figure 1: Segmented trailing edge flaps designations (top view). 


2.3. Wind tunnel data 

Two similar aircraft were wind tunnel tested in NASA Langley’s 12-foot Low Speed Tunnel. These aircraft allowed for 
the testing of several configurations, such that a baseline reference of data could be obtained. Between the two aircraft, 
a total of 275 static wind tunnel test runs were performed. 

One aircraft allowed for testing of a wing with full-span ailerons; each aileron surface encompassed one entire half 
wingspan. This same aircraft could also be equipped with conventional ailerons and flaps; each aileron surface 
accounted for one quarter wingspan, as did each flap surface. 

A second aircraft configuration enabled testing of the segmented trailing edge flaps, with each aileron segment 
independently controlled. With respect to the remaining control surfaces, the two aircraft were virtually identical, with 
the exception of a slightly larger elevator control surface on one of the two. 

The wind tunnel experiments were conducted at an average tunnel dynamic pressure of 2.95 lb/ft 2 . This correlates to a 
flow velocity of approximately 30 knots at sea level density. The typical airspeeds of the aircraft in question range from 
about 30 knots to 70 knots. A target airspeed of 45 knots was chosen as a reference for model development and data 
collection. This airspeed translates to a required lift coefficient of approximately 0.45 at an angle of attack of about 2 
degrees. Thus, all influence coefficient data collected is referenced to a nominal 2 degree angle of attack data set. 

2.3.1. Data extraction algorithms 


Table 1 : Matrix of wind tunnel run identifiers. 
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The tremendous amount of data available in the form of wind tunnel runs led to the development of algorithms to 
automatically extract the data of interest. After determining the location of this data in the data set, a matrix of wind 
tunnel run indices was developed, which provided a map to the desired information. Functions were developed for the 
elevator, the rudder, the segmented flaps, and sideslip perturbations that would call the run matrix to obtain the wind 
tunnel run numbers. These functions extracted the data in question from the wind tunnel runs and formed the non-linear 
controls influence functions of interest. Table 1 shows the wind tunnel run matrix used to extract the desired 
information. 

Some of the data extracted using this run matrix includes the reference coefficients for the aircraft as they vary with 
angle of attack and a unique algorithm was developed to extract this information. These reference coefficients are used 
to generate the traditional lift, drag, and moment profile curves for the aircraft. The influence functions can then be 
referenced with respect to these nominal values. 


3. MOTIVATION 


3.1. Influence functions 

Flexibility and stiffness influence coefficients are typically associated with vibrations or structural dynamics. These 
coefficients are used to describe how a system behaves when subjected to loads. 3 This is analogous to the aerodynamic 
and control derivatives used in aircraft dynamics. Typically, these coefficients are obtained at a reference or 
equilibrium point of interest of the system in question. 

The large number of control flaps on this aircraft presents a challenge in determining how to model the influence of 
each surface. Previous research on an aircraft with a similar control flap arrangement indicated that linear superposition 
could be used to model the total influence on the aircraft dynamics. 4 The influence functions developed here are 
implemented, using linear superposition, to create the overall system dynamics. 

Consider a vector containing the aerodynamic coefficients for the aircraft at an un-accelerated attitude with un-deflected 
(centered) control surfaces, fixed airspeed, angle of attack and sideslip angle. 

a„=[C z C D C r C, C m Cj 

At this reference attitude, it can be assumed that any changes in the elements of a 0 are due purely to control inputs. The 
vector of control inputs for this aircraft is as follows: 
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and t, e, la, ra, and r subscripts represent the throttle, elevator, left aileron flaps, right ailerons flaps, and rudder, 
respectively. For notational purposes, denote the reference coefficients and the control vectors simply as: 
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where the indices k= 6 and p= 19 in this example. 



Perturb or deflect the control surfaces some arbitrary amount away from the un-deflected position and obtain the change 
in the coefficients. The change in the i ,h coefficient away from its reference values due to this arbitrary set of 
control deflections can be defined as follows: 


A a,=a, -a, (13) 
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Consider the contribution from each control surface independently and expand the expression for the change in the 
coefficient as a linear superposition of the changes due to the p control flaps using partial derivatives, ignoring higher 
order terms: 
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Equation (14) provides the method for determining the resulting aerodynamic coefficients due to the control inputs: 
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The control inputs are referenced with respect to their un-deflected positions, 5- =0, therefore equation (15) 
becomes: 
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From equation (16), the influence coefficients are defined as the change in the aerodynamic coefficient a, resulting 
from a deflection of the control surface 5 , : 
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Equation (17) gives rise to some of the control derivatives traditionally used to describe aircraft dynamics. For 
example, the change in pitching moment coefficient due to a change in elevator control input is: 
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This specific control derivative is known as the elevator control power. 


Further, from equations (16) and (17), the linear controls influence functions can be defined as: 



(18) 



Traditionally, aerodynamic and control derivatives are collected for several flight conditions from either aircraft wind 
tunnel models, or from flight testing. Linearity is assumed and constant derivatives are used in describing the aircraft 
dynamics at or near the reference flight condition at which the derivatives were collected. This assumption can be 
considered valid only for small deflections of the control surfaces with small changes away from the reference flight 
condition. 


Aircraft control surfaces may exhibit non-linear influences on the aircraft dynamics. This can be due to asymmetries in 
the positive and negative deflections of a specific surface or from other phenomenon, such as changes in the flow due to 
the control deflection. In an effort to help make the controls influences more accurate, the wind tunnel data here was 
used to formulate non-linear influence functions for each control surface. These models are implemented in the aircraft 
dynamics simulation as functions that are called at every instant in the calculations to determine the appropriate control 
derivatives for the specified control deflections. 


The goal of this study was to develop functions that more accurately capture asymmetries in control deflections. Wind 
tunnel data was available that provided the capabilities to develop these functions for this specific aircraft. Trends were 
observed in the data and were used to develop function models that best represented each control input. The non-linear 
controls influence functions used to replace the linear functions are defined as: 
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where bjj are polynomials of degree n that represent best- fit trends to the wind tunnel data. Substituting equation (19) 
into equation (16) leads to the equation for the aerodynamic coefficients in terms of the non-linear controls influence 
functions. 
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3.1.1. Influence function models 

The first bits of information extracted from the available data were the aerodynamic force and moment coefficients for 
the aircraft. These coefficients describe the effects of angle of attack changes on the aircraft’s aerodynamic forces and 
moments/ The aerodynamic force coefficients are the lift, drag, and side-force coefficients seen in Figure 2. The 
moment coefficients are the rolling, pitching and yawing coefficients seen in Figure 3. The side-force, rolling moment, 
and yawing moment coefficients are small and have little influence on the aircraft dynamics with angle of attack 
changes, as can be seen in these figures. As a result, these coefficients are assigned zero reference values. The 
remaining coefficients are fitted using polynomials of the form: 


ci /((%) — d n cc H — + d^a -\-d\Oc-\~d 0 

The values of the coefficients for these polynomial fits are listed in Table 2. 
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Table 2: Reference aerodynamic coefficient polynomial fits. 


Aero. Coefficient 

d 6 

d 5 

d 4 

d 3 

d 2 

di 

d 0 

R 2 

Lift 

- 

- 

- 

-1.12E-04 

2.03E-04 

7.90E-02 

-3.64E-01 

0.999 

Drag 

- 

- 

- 

- 

9.48E-04 

-2.46E-05 

3.67E-02 

0.988 

Pitching Mom. 

-3.82E-07 

0.00E+00 

9.44E-06 

-3.25E-05 

-8.05E-04 

9.17E-04 

2.39E-02 

0.998 



Figure 2: Force coefficients versus angle of attack. Figure 3: Moment coefficients versus angle of attack. 


Changes away from the reference force and moment coefficients are due to the effects of the controls influence 
functions. Over 110 influence functions were extracted from the wind tunnel data and are used in describing the aircraft 
dynamics. Due to the large amount of data available and the limited space, only some of these representative influence 
functions are presented here. 

The controls influence functions, as described in equation (19), are modeled with polynomial fits of the form: 

bjj (Sj ) = d n S" + • • • + d 2 Sj + d ! Sj (22) 

The influence function models that describe how the elevator inputs affect the reference aerodynamic coefficients are 
presented in Table 3. 


Figure 4 and Figure 5 show the data used to fit these models. It is clear from the plots that three of the coefficients are not 
very sensitive to elevator deflections. These are the side-force coefficient, the rolling moment coefficient and the 
yawing moment coefficient. As a result, these are left as zero value influence functions. 



Figure 4: Elevator influence on the reference force coefficients. Figure 5: Elevator influence on the reference moment 

coefficients. 






Table 3: Elevator influence function models. 


Influence Function 
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- 
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- 
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- 
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0.00E+00 

- 


The models for the rudder effectiveness use the same influence function form of equation (22). These models are listed 
in Table 4. 


Table 4: Rudder influence function models. 


Influence Function 

d 3 

d2 

d, 

R 2 

bi,i9 

3.14E-07 

-8.02E-06 

-3.21 E-05 

0.379 

b2,19 
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5.85E-06 

8.83E-06 

0.913 

b 3 ,i9 

-9.03E-08 

-1.03E-05 

1.54E-03 
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1.02E-04 
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b5,19 

2.57E-07 

5.34E-06 

-2.51 E-05 

0.602 

b6,19 

- 

5.15E-06 

-6.60E-04 

0.999 


Two of the rudder controls influence functions models do not seem to fit the data as close as desired, although they do 
not produce significant changes in the coefficients. These two are b x 19 and b 5 19 , which correspond to changes in the 

lift coefficient and changes in the pitching moment coefficient, respectively, due to rudder deflections. Changes in the 
lift coefficient may be due to disturbances of the flow over the horizontal tail, possible resulting in a small loss of the 
total aircraft lift. Also, changes in pitching moment coefficient would likely be due to increased drag resulting from a 
positive or negative rudder deflection. These factors and aero-elastic effects may be provoking these erratic influences. 

Figure 6 and Figure 7 represent the data collected for the rudder influences. Evidence of asymmetric rudder influence 
can be seen in the function model for the yaw moment coefficient in Figure 7. Equal positive and negative deflections 
do not result in the same magnitude yawing moment coefficient. This asymmetry may be directly related to the 
asymmetric rudder geometry. 
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Rudder Influence on Moment Coefficients 


P -0.01 
-0.015 
- 0.02 


-20 


□ 

b 4,19 

❖ 

b 5,19 

0 

b 6,19 


4-- 

8 

9 

- 1 

— 

- i 

O _ 

- i 

9 


- i 


-10 0 10 
5 f (degrees) 


20 


Figure 6: Rudder influence on the reference force coefficients. 


Figure 7: Rudder influence on the reference moment 
coefficients. 





The segmented aileron flaps influence functions presented next are for flap number 15, the left wing outboard flap, flap 
number 11, one of the left wing mid flaps, and flap number 1, the left wing flap at the wing root. These three flaps 
exhibit the same trends in the data and can be used to show inconsistencies in the flow properties related to the 
aerodynamic coefficient influences. The influence function models for these flap segments are presented in Table 5. 


Table 5: Left wing sample flap (15, 11, and 1) influence function models. 


Influence Function 
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0.685 

b 2,10 

- 

6.38E-06 

4.68E-05 

0.625 

b3,10 

-1.15E-07 

6.35E-06 

-3.68E-05 

0.797 

b4,10 

-5.37E-08 

1.92E-06 

2.65E-04 

0.988 

bs.io 

-2.04E-07 

1.38E-06 

-2.36E-04 

0.941 

b&io 

- 

-1.18E-06 

-1.43E-05 

0.929 

bi,8 

3.63E-07 

9.49E-06 

7.10E-04 

0.956 

b2,8 

- 

3.67E-06 

2.30E-05 

0.689 

b3,8 

-1.82E-07 

-1.87E-06 

5.12E-05 

0.448 

b4,8 

-1.12E-08 

1.36E-06 

2.31 E-04 

0.991 

b5,8 

1.60E-07 

3.74E-06 

-3.43E-04 

0.940 

bd,8 

- 

-1.04E-06 

3.10E-06 

0.819 

bi,3 

5.01 E-07 

1.12E-05 

5.69E-04 

0.882 

b2,3 

- 

3.10E-06 

1.45E-05 

0.477 

b3,3 

-3.1 8 E-07 

1.2 IE-06 

4.27E-05 

0.911 

b4,3 

1.1 IE-07 

1.35E-06 

7.06E-06 

0.937 

b5,3 

1.39E-08 

8.54E-07 

2.74E-04 
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It is clear from the data presented in Figure 8 - Figure 13 that linear fits to this data may not accurately describe the 
influence of these flaps on the dynamics of the aircraft. The segmented flaps affect the lift coefficient in a general cubic 
curve-fit manner; the influences dropping off with increased deflections from nominal zero deflection. The contribution 
of these flaps to the drag coefficient shows a quadratic trend, where both positive and negative deflections provoke an 
increase in drag. The side-force coefficient, as expected, is not very sensitive to the deflection of the segmented aileron 
flaps. 
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Figure 8: Flap number 15 effect on the reference force 
coefficients. 
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Figure 9: Flap number 15 effect on the reference moment 
coefficients. 




Figure 10: Flap number 1 1 effect on the reference force 
coefficients. 



Figure 12: Flap number 1 effect on the reference force 
coefficients. 



Figure 1 1 : Flap number 1 1 effect on the reference moment 
coefficients. 



Figure 13: Flap number 1 effect on the reference moment 
coefficients. 


The moment coefficient influence functions for these flaps show intuitive behavior and some possible flow 
inconsistencies. Positive deflections of these flaps should result in an increase in the rolling moment coefficient, while 
negative deflections should result in a decrease in this coefficient. The opposite is true for the right wing flaps. This 
trend is seen in the data presented in Figure 9, Figure 11, and Figure 13. It is interesting to note that this influence on the 
rolling moment shows another expected trend where the outboard flap, flap number 15, has a greater impact on the 
rolling moment than the inboard flap, flap number 1. Previous research, involving an aircraft with similar control 
surfaces, indicated this same trend with the inboard control flaps having relatively low impact on the roll moment. 6 

The yawing moment coefficient is not particularly sensitive to these flaps deflections but does show an expected trend. 
Deflections of these flaps, positive or negative, serve to increase the drag produced by that wing, causing an increase or 
decrease in the yawing moment. From this data, deflection of the left wing flaps will generally cause a decrease in the 
yawing moment coefficient, which is expected. 

The pitching moment coefficient influence functions show some unexpected behavior that is consistent between the left 
and right wing segmented flaps. Figure 9 shows an increase in the pitching moment coefficient due to a negative 




deflection of flap segment number 15, and vice versa for a positive flap deflection. This is what is generally expected 
since a negative deflection of the flaps is a trailing edge up deflection and that would typically cause an increase in the 
pitching moment, or a nose-up pitching of the aircraft. As we look at flap number 1 1, moving along the left wing from 
the wingtip to the wing root, we see the same behavior or trend. Comparing this behavior to Figure 13, the influence of 
flap segment number 1, we see a reversal in the pitching moment influence. This reversal in behavior is evident in the 
two inboard flap segments for each wing. Since these are the closest flaps to the fuselage, flow interactions may be 
occurring that result in this behavior. 

The flap segment influence functions not presented here showed similar behavior to flaps 15 and 11. All influence 
functions for the segmented trailing edge flaps have been implemented in the aircraft dynamics simulator. 

Influence functions were also created for the sideslip effects on the force and moment coefficients. These influence 
function models were simple models because data was only available for -4 degrees, zero degrees, and 4 degrees 
sideslip deflections. 

At this stage in the simulation development, influences on the aerodynamic coefficients due to throttle inputs will be 
assumed zero. This is not to say that these would not have an impact on the reference coefficients. In fact, this aircraft 
is a puller propeller aircraft and changes in throttle would lead to possible increases in lift coefficient due to additional 
flow velocity over the aircraft. Also, engine torque would likely influence the rolling moment coefficient. 


3.1.2. Additional models for simulation 

The influence functions developed here are static functions that describe the aircraft static stability. Additional 
modeling was required to ensure an accurate representation of the aircraft dynamics. In particular, damping was 
required in areas that did not yet have damping. Additionally, an engine model was needed for the aircraft. 

Table 6 is a list of the stability derivatives or influence coefficients typically used to model aircraft dynamics. Some of 
these derivatives are already embedded in the model aircraft in the form of the influence functions. The lift-curve slope, 
for example is taken into account in the models extracted from the wind tunnel data for the reference aerodynamic 
coefficients. 


Table 6: Longitudinal and lateral/directional aircraft stability derivatives. 



Longitudinal Stability Derivatives 


Lateral/Directional Stability Derivatives 

Derivative 

Description 

Modeled As 

Derivative 

Description 

Modeled As 


lift-curve slope 

ref. coefficient models 

c '„ 

dihedral derivative 

influence function 

c„ a 

pitch stiffness 

ref. coefficient models 

C n> 

yaw stiffness 

influence function 

C„ q 

pitch damping 

estimation model 

C, p 

roll damping 

estimation model 

c my 

tuck derivative 

-0 (subsonic) 

c„. 

yaw damping 

estimation model 

C m , 

alpha-dot derivative 

-0 (alpha-dot small) 

c „ p 

yaw moment due to roll rate 

~0 (assumed small) 

C Dr 

speed damping 

-0 (low speed) 

C,, 

roll moment due to yaw rate 

-0 (small yaw rate) 

c °„ 

drag-curve slope 

ref. coefficient models 

Cj, 

sideforce due to sideslip 

influence function 

C w 

lift vs. speed slope 

-0 (low speed) 


sideforce due to yaw rate 

-0 (small yaw rate) 

c La 

acceleration derivative for lift 

~0 (alpha-dot small) 

Cj, 

sideforce due to roll rate 

-0 (assumed small) 

C L, 

pitch-rate dependent lift 

~0 (q small) 

c 

yaw mom. due to sideslip rate 

~0 (small sideslip 
rate) 




C 

sideforce due to sideslip rate 

~0 (small sideslip 
rate) 





roll moment due to sideslip 

-0 (small sideslip 




c ‘i 

rate 

rate) 


Many of the derivatives in Table 6 apply to high speed or high maneuverability aircraft and can be assumed to be zero. 
The aircraft in question is a low speed aircraft and, for the purposes of developing this simulation, it was assumed that 






the maneuvers performed are low speed or fairly benign. This assumption eliminates many of the derivatives from the 
list. 


Given the assumptions of this model, three additional stability derivatives were modeled. These are the roll damping, 
pitch damping, and yaw damping derivatives. The pitch damping and yaw damping derivatives are statistical models 
that depend on the geometric properties of the aircraft. 1 ' 7 The roll damping derivative is an approximated guess that 
falls within historical levels and produces realistic dynamic behavior when compared to observed aircraft dynamics. 
These three damping derivative values are listed in Table 7. 


Table 7: Damping derivatives. 


Derivative 

Description 

Est. Value 

c, p 

roll damping 

-0.25 

C mq 

pitch damping 

-5.31 

c„ 

yaw damping 

-0.50 


One additional model was developed for use in simulating the aircraft dynamics. Engine performance data was not 
available for this aircraft and a simple model was developed for the internal combustion, propeller driven engine. An 
estimate was made on the maximum static thrust of the engine and an exponential decay model was used to simulate a 
decrease in the available thrust with increasing airspeed. 

In general, an aircraft dynamics simulation would also require an altitude density model. However, for purposes of the 
research at hand, all flights can be assumed to take place at or near sea level conditions, thus eliminating the need for 
this model. 


4. RESULTS 


4.1. Simulation Capabilities 

The aircraft dynamics simulation developed using these non-linear controls influence functions results in a multiple 
capability tool. Of particular interest and usefulness is the ability to understand the dynamics of an aircraft that 
implements innovative controls surfaces for multifunctional purposes. Computer simulations allow for scenarios to be 
run and controllers to be tested prior to actual flight testing. 

The simulation developed here will serve as a basis for continuing research. This model can easily be modified to 
account for changes in aircraft design by inputting new data obtained from either wind tunnel experiments or CFD 
analysis. Also, this model can easily be modified to incorporate aero-elastic effects. 




Figure 14: Dynamics visualization window. 


Figure 15: Aircraft trajectory visualization window. 





Another great capability of this aircraft dynamics model is the ability to extract linear plant models. A reference flight 
condition of interest can be chosen and a canonical form system model extracted for use in controller development. 
These linear models allow for the determination of system stability and frequency domain analysis. Control law 
developed using these linear models can easily be tested on the full non-linear aircraft model. 

Free or controlled responses can be visualized via a 3-dimensional aircraft model translating and rotating in 3- 
dimensional space as seen in Figure 14. The trajectory traveled by the aircraft during a specified maneuver can be seen 
in a similar 3 -dimensional manner, as in Figure 15. Finally, the aircraft state time histories resulting from a maneuver 
can be analyzed from plots as can be seen in Figure 16. 



Figure 16: Aircraft states time histories (seconds). 


In understanding the dynamics of an aircraft it is important to quantify its performance capabilities. The simulation 
developed here allows for the generation of a performance envelope for the aircraft, given the assumptions made in 
characterizing the dynamics. Information such as stall speed, control surface trim conditions, climb rate, and more can 
be extracted from data obtained from the simulation. 

5. CONCLUSIONS 

Non-linear controls influence functions have been implemented in a simulation to characterize the dynamics of an 
aircraft with novel controls effectors. The segmented trailing edge flaps on the aircraft in question provide additional 
degrees of freedom in control that could not be captured by traditional aircraft modeling techniques. The influence 
function method implemented in this study allows for control input asymmetries to be accurately represented in the 
aircraft dynamics. 


As a result of this study, a dynamics model has been developed that allows for extraction of linear plant models and 
more. The model can be adapted with new influence functions if new wind tunnel data becomes available for the same 
aircraft or for an aircraft with a similar controls arrangement. The free or controlled response of the aircraft non-linear 




dynamics can be visualized in a 3 -dimensional environment, in addition to presentation of the aircraft trajectory and the 
aircraft states time histories. Additional possibilities of this dynamics simulation can be endless and are left to the 
imagination of the reader. 
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